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ABSTRACT 

Transit timing analysis may be an effective method of discovering additional bodies 
in extrasolar systems which harbour transiting exoplanets. The deviations from the 
Keplerian motion, caused by mutual gravitational interactions between planets, are 
expected to generate transit timing variations of transiting exoplanets. In 2009 we 
collected 9 light curves of 8 transits of the exoplanet WASP-lOb. Combining these data 
with published ones, we found that transit timing cannot be explained by a constant 
period but by a periodic variation. Simplified three-body models which reproduce the 
observed variations of timing residuals were identified by numerical simulations. We 
found that the configuration with an additional planet of mass of ~0.1 Mj and orbital 
period of ^5.23 d, located close to the outer 5:3 mean motion resonance, is the most 
likely scenario. If the second planet is a transiter, the estimated flux drop will be ^0.3 
per cent and can be observable with a ground-based telescope. Moreover, we present 
evidence that the spots on the stellar surface and rotation of the star affect the radial 
velocity curve giving rise to spurious eccentricity of the orbit of the first planet. We 
argue that the orbit of WASP-lOb is essentially circular. Using the gyrochronology 
method, the host star was found to be 270 ± 80 Myr old. This young age can explain 
the large radius reported for WASP-lOb. 

Key words: planetary systems - stars: individual: WASP-10. 



1 INTRODUCTION 

Analysing transit timing variations (TTV) of exoplanets is 
expected to be an efficient method to d i scover add i tional 
low-mass planets l|Miralda-Escudel |2002| ; ISchneiderl 120041: 



* Partly based on observations made with the 0.6 and 2.0-m tele- 
scopes of the Rozhen National Astronomical Observatory, which 
is operated by the Institute of Astronomy. Bulgarian Academy 
of Sciences and the 90-cm telescope of the University Observa- 
tory Jena, which is operated by the Astrophysical Institute of the 
Friedrich Schiller University, 
f E-mail: gm@astro.uni-jena.de 



Holman fc Murravl 120051 : lAgol et al.l 120051 : ISteffen et all 
20071 ). Many studies have been performed to detect TTV 
signals but they resulted only in constr aints on the parame - 
ters of hy pothetical second planets fe.g. lGibson et al]|2009l ). 
Recently, Maciciewski et al. (2010) detected variation in the 
transit timing of WASP-3b. They found that a configuration 
with a hypothetical second planet of mass of ~15 Mm, lo- 
cated close to an outer 2:1 mean motion resonance (MMR) 
may reproduce the observed TTV signal. iLendl et al] (|2010h 
found some indications for the TTV for OGLE2-TR-L9b but 
no preliminary solution was proposed. 

WASP-lOb is an exoplanet discovered by 
IChristian et~ai1 (120091 ). The mass and radius of the 



2 G. Maciejewski et al. 



Table 1. List of telescopes: FoV is the field of view of the instrument and Ntr is the number of observed transits. 



ID 


Telescope 


Observatory 


Detector 


A?tr 






Location 


CCD size 


Field of view 


1 


0.6-m Cassegrain 


Astronomical Observatory, N.Copernicus University 


SBIG STL-1001 


1 






Piwnice near Torun, Poland 


1024 X 1024, 24 /im 


ll.'S x 11.'8 


2 


0.6-m Schmidt* 


University Observatory Jena, 


CCD-imager STK 


2 






GroBschwabhausen near Jena, Germany 


2048 x 2048, 13.5 (im 


52.' 8 x 52f8 


3 


0.6-m Cassegrain 


National Astronomical Observatory 


FLI PL09000 


4 






Rozhen, Bulgaria 


3056 x 3056, 12/itn 


17.'3 x 17.'3 


■1 


2-m Ritchey-Chretien 


National Astronomical Observatory 


PI Vers Array :1300B 


1 






Rozhen, Bulgaria 


1340 x 1300, 20/im 


5f8 x 5.'6 


5 


0.8-m Tenagra II 


Tenagra Observatories 


SITe 


1 






Arizona, USA 


1024 X 1024, 24/^m 


14.'8 x 14.'8 


6 


1.5-m Ritchey-Chretien 


Gunma Astronomical Observatory 


Andor DW432 


0** 






Takayama, Japan 


1250 X 1152, 24^m 


12.'5 x 11.'5 



no scientific output due to bad weather conditions 



planet were found to be 2.mt° ^ Mj and 1.28±g;g| Rj, 
respectively, larger than interior models of irradiated giant 
planets predict. The planet orbits its host star in ~3.09 
d at an orbital semi-major axis of 0.03691q'qqj4 au. The 
eccentricity of the orbit was reported to be 0.059to:oo4 
- a l arge value for a clos e-in planet. In the discovery 
paper IChristian et~all (|2009r i determined transit depth and 
duration to be 29 mmag and 2.36 h, respectively. 



IJohnson et all (|2009T l and Ijohnson et all (|201(t) rede- 
termined the planet's parameters with a high-precision light 
curve using the University of Hawaii 2. 2-m telescope with 
a photometric precision of 4.7 x 1CP 4 and the mid-transit 
time error of 7 s. The mass of the planet was found to be 
3-151q'Jj Mj and the planetary radius turned out to be no- 
ticeably smaller, i.e. 1.08=1=0.0 2 Rj. The latter result was not 
confir me d by the studies of iKreicova, Budai fc Krushevskal 
l|20ldl or iDittmann et alj i|2010T ). The first team found the 
radius of WASP- lOb to be in exc e llent agreement with the 
value obtained bv lChristian et ab I il2009T ), and a similar con- 
clusion was reached by Dittmann et al.l (|2010h with no TTV 
signature detected. 

The host star (V = 12.7 mag) is a K5 dwarf with effec- 
tive t emperature of 4675±1 K, located 90±20 pc from the 
Sun l|Christian et al l I2009T I. ISmith et al l l|2009T ) confirmed 
previous evidences that the light curve of WASP- 10 exhibits 
a rotational photometric variability caused by star spots. 
The stellar rotation period was foun d to be 11.91 ± 0.05 d - 
in agreement with ~12 d reported bv lChristian et al l (2009). 

The presence of an additional small planet was 
suggested to p ump the eccentricity of WASP- 10b 
(IChristian et al. I I2009T ). This system was therefore a 
good target for TTV studies. 

The order of this paper is as follows. In Section 2 obser- 
vations and data reduction are summarised. The discovery of 
the periodic transit timing variation is presented in Section 
3. The process of identifying and verifying of the most likely 
2- planetary model of the WASP- 10 system is described in 
Section 4. A discussion of the proposed scenario is presented 
in Section 5. Finally, conclusions are collected in Section 6. 



2 OBSERVATIONS AND DATA REDUCTION 

We collected 9 light curves of 8 transits of WASP-lOb dur- 
ing a dedicated international observing campaign involving 
telescopes spread worldwide at different longitudes (listed 
in Table [T] and 0). The transit on July 31, 2009 was ob- 
served with two instruments in different observatories. The 
telescope diameters of 0.6 to 2.0 m allowed us to collect pho- 
tometric data with 1.1 — 2.7 mmag precision. Observations 
generally started ~1 hour before the expected beginning of a 
transit and ended ~1 hour after the event. However, weather 
conditions and schedule constraints meant that we were not 
always able to follow this scheme. 

Standard IDL procedures (adapted from DAOPHOT) 
were used for the reduction of the photometric data col- 
lected at the National Astronomical Observatory at Rozhen, 
Bulgaria, and computing the diffe rential ape r ture p hotome- 
try. Using the method of I Everett fc Howelll (|200lT ) several 
stars (4 to 6) with photometric precision better than 5 
mmag were selected to create an artificial standard star 
used for differential photometry. Data from the remaining 
telescopes were reduced with the software pipeline devel- 
oped for the Semi-Automatic Variabilit y Search sky survey 
(Nicdziclski, Macicicwski & Czart 2003). To generate an ar- 
tificial comparison star, 30-50 per cent of the stars with the 
lowest light-curve scatter were selected iteratively from the 
field stars brighter than 3 mag below the saturation level. 
To measure instrumental magnitudes, various aperture radii 
were used. The aperture which was found to produce light 
curves with the smallest scatter was used to generate a final 
light curve. 



3 RESULTS 

3.1 Light curve analysis 

A model-fittin g algorithm available via the Ex oplanet Tran- 
sit Database (|Poddanv. Brat fc Peichal l201fj ) was used to 
derive transit parameters: duration, depth and mid-transit 
time. The procedure e mploys the OCCULTSMALL routine of 
iMandel fc Agoll |2002l ) and the Levenberg-Marquardt non- 
linear least squares fitting algorithm which also provides 
uncertainties. As our research was focused on determining 
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Figure 1. Light curves of WASP-lOb transits in individual runs. The best-fitting models are shown as continuous lines. The transit on 
31 July was observed simultaneously by two observatories. 
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Figure 2. The observation minus calculation (O — C) diagram for W ASP- 10b generated f or linear ephem e ris ba sed o n literature data 
(continuous line). Open circles denote literature data taken from [ch ristian et al.l j200Sl) ■ Ijohnson et al.l ll2010f ) and iDittmann et al] 
Il2010h . Open squares are based on reanalysed photometry from lKreicova eT'aTTilijOlOl '). Filled symbols denote results from our campaign. 
A significantly be tter fit may be obtaine d for the new ephemeris given by Eq. 1 (dashed line). 3cr error bars were taken for mid-transit 
time reported by IDittmann et all 1120101) . Their light curve was not corrected for a linear trend which is well visible in out-of-transit 
phases (J. Dittmann 2010, priv. comm.). Our tests showed that this effect significantly affected the accuracy of the mid-transit time. 
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Table 2. The summary of observing runs: ID - the identification 
of the instrument according to Table [T] jV ex p - number of useful 
exposures, T exp — exposure times. Dates UT are given for the 
beginning of nights. 



Run 


Date 


ID 


Filter 






1 


2009 Jan. 5 


1 


R 


254 


50 


2 


2009 July 31 


2 


R 


144 


60, 80 


3 


2009 July 31 


3 


R 


116 


120 


4 


2009 Aug. 3 


3 


R 


122 


90 


5 


2009 Aug. 28 


3 


R 


99 


90 


6 


2009 Aug. 31 


2 


R 


203 


50, 60 


7 


2009 Sep. 3 


3 


R 


105 


90 


8 


2009 Oct. 1 


4 


V 


548 


15 


9 


2009 Nov. 17 


5 


R 


76 


90 



mid-transit times, the number of parameters to be fitted 
could be reduced. An impact parameter b = acosi/P* = 
0.299^0 043, where a is semi-major axis, i is inclination 
and R , is host-star radius, was taken from Ijohnson et al.l 
(2009) and was fixed during the f itting procedure. W e used 
the linear limb-darkening law by IVan Hammel (|l993h with 
the limb-darkening coefficient linearly interpolated for the 
host star in a given filter. A first- or second-order poly- 
nomial was used to remove trends in magnitude or flux. 
The mid-transit times were corrected from UTC to Terres- 
trial Dynamical Time (TT) and t hen transformed into BJD 
ijEastman. Siverd fc Gaudl l2010t l. Light curves with best- 
fitting models and residuals are shown in Fig. [1] and derived 
para meters are collected i n Table [3] 

iKreicova et alj (|2010t ) also published observations and 
P-band light curves of four complete transit events of 
WASP-lOb. In their analysis they focussed mainly on the 
determination of the planet radius and did not study the 
mid-transit times, O — C diagram or TTV. We determine 
4 mid-transit times for these observations and placed them 
into the context of our data to create a larger sample for 
TTV analysis. We used the same method to determine the 
mid-transit times from their light curves as for all other ob- 
servations presented here. Corresponding, mid-transit times 
are presented in Table [3] 

Across all data sets, two transits on epochs 222 and 232 
were observed simultaneously by two different telescopes. 
The data were reduced with different pipelines and by dif- 
ferent teams. The transit on epoch 222 was observed with 
60-cm telescopes in Jena and at Rozhen. The difference in 
mid-transit times was found to be ~11 s. The transit on 
epoch 232 was covered by observations in Jena and the 
50-cm telescope in Stara Lesna in Slovakia (|Kreicova et al.l 
2010). In this case the discrepancy of timing was found to 
be ~18 s. Timing differences are well within error bars of 
individual determinations for both transits. This practical 
test strengthens the reliability of our mid-transit time de- 
terminations. 



3.2 Detection of a TTV 

While determining new linear ephemeris, we noted that a 
linear fit of the epoch E and orbital period Pb of WASP- 
lOb resulted in reduced Xrcd = 13.8. Individual mid-transit 
errors were taken as weights in the fitting procedure. Such 



a high value of Xred suggests the existence of additional sig- 
nal in the O — C diagram (Fig. [2l). The Lomb-Scargle peri- 
odogram (|Lombll976J : IScarglelll982l ) of the residuals revealed 
the existence of 2 peaks of similar significance and frequen- 
cies of 0.175 and 0.183 cycl P" 1 (Fig.©. For each frequency, 
an ephemeris was refitted with the linear trend plus a sinu- 
soidal variation in the form: 

T b = To + E x P b + Auv sin ^2tt - ~ t f ttv ) (1) 

where To is the mid-transit time for initial epoch (E = 0), 
Attv is the semi-amplitude of detected transit timing varia- 
tion, Ettv is the epoch offset of the TTV signal and P ttv is 
its period. The minimal x? e d = 3.2 was obtained for /* tv = 
0.183 cycl P" 1 and resulted in T = 2454357.86011±0.00047 
BJDtt, P b = 3.0927183 ± 0.0000021 d, A ttv = 0.00120 ± 
0.00036 d, P ttv = 1.8 ± 0.4 and Pf v = 5.473 ± 0.012 P b . 

The procedure run for /| tv = 0.175 cycl P^ 1 gave 
X? cd = 4.1, P 2 ttv = 5.7172 ± 0.0082 P b and similar values 
of To and Pb. In this case the fit turned out to be poorer, 
thus we used the ephemeris based on parameters obtained 
for P? v . 

The ephemeris calculated according to Eq. 1 (and P" v ) 
is plotted in Fig. [5] with a dashed line. The O — C values 
which are used in further analysis were calculated according 
to the linear part of Eq. 1. 

To check if a random distribution of data points in the 
O — C diagram may favour the detected frequencies, 10 
fake data sets were generated. The original residuals were 
replaced with random values with a white-noise distribu- 
tion of the amplitude. Then, the Lomb-Scargle algorithm 
was run to find a dominant frequency in each fake diagram. 
The histogram of these frequencies, with the number of bins 
equal to the square root of the number of data points in the 
sample, is presented in Fig. [4] Detected TTV frequencies, 
/" v and /2 tv , are not associated with any significant peak 
in t he histogram. 

iKipping fc Bakosl (|2010al ) pointed out that the transit 
phasing which is the time difference between the expected 
mid-transit moment and the nearest data point in a light 
curve may generate a spurious TTV signal. For example, 
for a light cu rve of 60-s cadence time difference is expected 
to be ±30 s l|Kipping fc Bakosll2010bh . Cadences of our in- 
dividual light curves are in a wide range between 19 and 
143 s (see Table [3} with the median value of 63 s. This 
value is ~3 times smaller than the amplitude of the detected 
TTV signal. As it is shown in Fig. [5] no correlation between 
individual residuals and cadences was detected. A formal 
least-square fit to this data set resulted in the correlation 
coefficient R = —0.36 which clearly reveals no correlation 
at the relevant level. The checks discussed above strengthen 
the detection of the TTV signal in our data set. 



4 TWO-PLANETARY MODEL 

The possible non-z ero eccentricity postulated by 
IChristian et al.1 l|2009h might potentially have resulted 
from confusion with a two-planet system in which both 
planets or bit their host star in circular o rbits in an inner 2:1 
resonance l|Anglada-Escude et alj|2010h . Although we argue 
later that eccentricity of WASP-lOb is indistinguishable 
from zero, the case of an inner perturber was considered 
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Table 3. Parameters of transit light-curve modelling. To denotes the mid-transit time, Tj - transit time duration, 8 
- depth, a - averaged standard deviation of the fit and E — epoch. O — C values were calculated according to the 
linear part of Eq. [l] and the best-fitting parameters. O — C values are given both in days and in units of multiplies 
of the uncertainty of the mid-transit times. t cad is the mean cadence of a light curve. Results of reanalysed data from 
iKreicova et al.l J2010l> are collected below the horizontal line. BJD times are based on Terrestrial Dynamical Time (TT). 



Run 


To 


T d 


6 


<7 


E 


O-C 


O-C 


^cad 




BJD TT - 2450000 


(min) 


(mmag) 


(mmag) 




(d) 


(To errors) 


(s) 


1 


4837.23116 ± 0.00037 


128.9 ±1.4 


34.1 ± 0.8 


2.7 


155 


-0.00030 


-0.8 


56 


2 


5044.44427 ± 0.00032 


128.9 ±1.2 


33.5 ± 1.1 


1.7 


222 


±0.00069 


+2.2 


87 


3 


5044.44415 ± 0.00040 


129.1 ± 1.5 


33.1 ± 0.7 


2.6 


222 


±0.00057 


+ 1.4 


130 


4 


5047.53696 ± 0.00028 


129.4 ± 1.0 


32.4 ±0.5 


2.0 


223 


±0.00066 


+2.3 


103 


5 


5072.27861 ± 0.00053 


131.1 ±2.1 


35.2 ± 0.8 


2.5 


231 


±0.00056 


+ 1.1 


102 


6 


5075.37129 ± 0.00024 


129.2 ±0.8 


33.0 ±0.7 


2.0 


232 


±0.00053 


+2.2 


63 


7 


5078.46485 ± 0.00044 


130.4 ±1.5 


34.5 ± 1.1 


2.5 


233 


±0.00137 


+3.1 


102 


8 


5106.29681 ± 0.00021 


129.8 ±0.7 


34.4 ±0.7 


1.1 


242 


-0.00115 


-5.4 


19 


9 


5152.68923 ± 0.00039 


128.8 ± 1.3 


33.3 ±2.4 


2.1 


257 


±0.00051 


+ 1.3 


143 




4775.37834 ± 0.00042 








135 


±0.00125 


+2.9 


35 




5075.37150 ± 0.00034 








232 


±0.00073 


+2.2 


26 




5109.39135 ± 0.00052 








243 


±0.00068 


+ 1.3 


36 




5112.48528 ± 0.00035 








244 


±0.00189 


+5.4 


36 
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Figure 3. The Lomb-Scargle periodogram generated for timing 
residuals plotted in Fig. [2] showing the existence of a periodic 
signal. The most significant peaks are indicated. The upper limit 
of the periodogram is determined by the Nyquist frequency for 
the data set. The vertical dashed line marks the rotational period 
of the host star (see Section 5). 

for completeness of our analysif]. Assuming circular orbits, 
the mass of the perturbing planet M c depends on the 
mass of the outer, more mas sive transiting planet Mb an d 
its apparent eccentricity eb (|Anglada-Escude et al-lfeoiOh . 
This results in M c « 0.14 Mj in the WASP-10 system. 
Such a planet is expected to produce strong gravitational 
perturbations which should be visible as the TTV of 
WASP-lOb. To check this scenario, we generated synthetic 

— C diagrams for WASP-lOb in systems with a second 
planet. Parameters o f the transiting plan et and its host 
star were taken from I Johnson et al.l (2009). Initial circular 
and coplanar orbits were assumed for both planets. The 
perturbing planet was put in an orbit with the semi-major 

1 We do not analyse the scenario with an exomoon orbiting 
WASP-lOb be cause no transit duration variation, predicted by 
iKippind (|2009h . was detected in our data. 
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Figure 4. The distribution of dominant frequencies for 10 5 fake 
O — C diagrams with observed data points replaced by the white 
noise. The number of bins is equal to the square root of the num- 
ber of data points in the sample and the width of individual bins 
is 0.001 cycl P" 1 . Detected TTV frequencies, /* tv and f* tv , are 



indicated. No significant peak is associated with any of them what 
indicates that a spectral window is not a source of detected peri- 
odicities. 



axis between 0.0218 and 0.0257 au (±0.0020 au away from 
the 2:1 orbital resonance). Calculations wer e performed us- 
ing the Mercury package (|Chambers!ll999l ) employing the 
Bulirsch-Stoer integrator. Simulations covered 270 periods 
of WASP-lOb, i.e. the time span covered by observations. 
No configuration was found to reproduce periodicity close 
to P* tv and P 2 ttv . The procedure was repeated for inner 
3:1 and 3:2 MMRs but these cases also brought negative 
results. 

To search for configurations with an outer perturber re- 
producing Pl tv and P" v > synthetic O—C diagrams were gen- 
erated with the ptmet code which is based on pertur bation 
theory i|Nesvornv fc Morbidellill2008l : H esvorny 2009). Both 
planets were put in initial circular orbits. The semi-major 
axis a c of the perturber varied between 0.0500 and 0.3000 au 
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Figure 5. Absolute values of individual residuals in transit timing 
vs. cadence. Symbols denote the same data sets as in Fig. [2] No 
correlation between both quantities was found. 
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Figure 7. The Lomb-Scargle periodogram for the radial velocity 
residuals after removing the transiting planet in a circular orbit. 
A peak corresponding to the stellar rotation is marked. 



Table 4. Outcr-pcrturbcr solutions which reproduce the ob- 
served O — C variation. p ttv indicates which periodicity in the 
O — C diagram (P" v or P^ tv ) is reproduced by a solution, a c 
denotes the semi-major axis of the perturbing planet, M c is its 
mass, P c is its orbital period, K c is the expected semi-amplitude 
of the radial-velocity variation and Xr e d i s the lowest value of 
reduced chi-square for direct model fitting. 



No. 


pttv 




M c 


Pc 


K c 


Y 2 
^red 






(au) 


{M 3 ) 


(d) 


(m s- 1 ) 


1 


1 


0.0536 


0.10 


5.2293 


14.2 


1.5 


2 


2 


0.0539 


0.10 


5.2647 


14.1 


2.5 


3 


2 


0.0682 


0.55 


7.4962 


69.1 


2.8 


1 


1 


0.0686 


0.55 


7.5677 


68.9 


2.8 



in steps of 0.0001 and 0.0002 au for a c ^ 0.1 and a c > 0.1 au, 
respectively. The amplitude of the TTV signal scales nearly 
linea rly with the perturber's mass |Nesvornv fc Morbidellil 
2008), thus M c was fixed and set equal to 0.3 Mj for the 
preliminary identification. Each simulation covered 270 pe- 
riods of WASP- 10b. Periodic signals reproducing P" v and 
Pi v were identified close to 5:3 and 5:2 MMRs for perturber 
masses of ~0.1 and ~0.5 Mj, respectively. 

To refine parameters of a perturbing body, the synthetic 
O — C diagrams were calculated with the Mercury code 
and the Bulirsch-Stoer integrator. The initial parameters 
of planet c were taken from previous simulations and then 
refined with the mass varying in steps of 0.05 Mj. Each sim- 
ulation covered 31000 days (i.e. ~10000 periods of WASP- 
10b). The synthetic O — C diagrams were directly fitted to 
observed data by shifting along the time axis and adjusting 
to the time span covered by observations. The best-fitting 
solutions are collected in Tabled The chi-square test favours 
solution 1 which reproduces Pi tv with minimal Xred = 1-5- 
This value is significantly smaller than Xred obtained in an 
analogous way for remaining configurations. The O — C dia- 
gram with solution 1 is presented in Fig. [6] where the resid- 
uals are also plotted. 



5 SPECTROSCOPIC REANALYSIS AND 
STELLAR ACTIVITY 

To re-assess the possible solution s, we reanalyse d radia l 
velocity (RV) data published by IChristian et al. 12009). 
We u sed the Systemic Console software (|Meschiari et al.l 
2009). The host star is known to be varia ble due to stellar ro- 
tation and starspots |Smith et alj|2009l ). It was shown that 
this effect may generate the RV variations up to a few hun- 
dred m s" 1 l|Hatzesll2002T) and mimic a p lanetary companion 
on a Keplerian orbit ( Desort et al.ll2007l ) . We reanalysed RV 
data assuming WASP- 10b has a circular orbit. This resulted 
in a best-fitting model with Xred = 3-9 and rms — 49.6 
m s _ . The periodogram of the residuals (Fig. [7| reveals a 
periodicity of ~ 11.84 d - a value close to the per iod of stellar 
rotation P r ot = 11.91 ± 0.05 d determined by ISmith et all 
(2009). The fitting procedure was repeated with an addi- 
tional sinusoid signal of period P ro t and a floating phase 
and amplitude. The eccentricity of WASP- 10b, eb, was also 
allowed to vary. This resulted in a significantly improved 
model with Xred = 2.5 and rms = 31.4ms -1 . The eccentric- 
ity was found to be eb = 0.013 ± 0.063 and it is statistically 
indistinguishable from zero. Therefore, we adopted eb = 0.0 
in further analysis. This result is not surprising as tidal in- 
teractions with the host star ar e expected to circularise the 
orbit of the planet (jZahnlll977l V 

To check if the perturbing planet may be hidden in the 
RV residuals of the model composed of the planet b and 
starspots, the second planet with physical and orbital pa- 
rameters as in solutions 1-4 was put into the system. A cir- 
cular Keplerian orbit of the second planet and coplanarity 
were assumed for simplicity. We noted that only solutions 1 
and 2 significantly decrease the rms of the RV model to 24.7 
and 22.5 m s _1 , respectively. Combining results of the TTV 
and RV analysis, we found solution 1 the most likely. The 
RV curves are plotted in Fig. [8] for individual components, 
i.e. WASP-lOb, WASP-lOc and the starspots effect. 

The RV semi-amplitude generated by stellar rota- 
tion was found to b e ~65 m s _1 . RV observations by 
IChristian et~ai1 |2009ri span from August 2007 to Febru- 
ary 2008 and are partially covered b y photometric obser - 
vations from the S uperWASP survey ( Pollacco et al.l l2006'). 
ISmith et all (2009) determined the amplitude of photomet- 
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Epoch 

Figure 6. The O — C diagram for WASP-lOb with the best-fitting model (solution 1 in Tablef!}. The mass of the perturber is 0.10 Mj 
and its orbital period is ~5.23 d. Symbols are the same as in Fig. [2] Residuals are plotted in the bottom panel. 
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Figure 8. Radial-velocity curves after removing remaining 2 com- 
ponents for a) WASP- lOb, b) WASP-l O c and c) starspots effect. 
Data were taken from IChristian et al. Two points with 

large error bars were obtained in poor weather conditions and 
ma y exhibit systematic errors due to scattered moonlight (see 
e.g. lLatham et alj|2010h . 



ric variation in the 2004 and 2006 seasons only. To find 
this variation in t h e 200 7 season and verify determinations 
from ISmith et al] l2009h. the analysis of variance method 
(ANOVA, ISchwarzenberg-Czernvl 1 1996h was applied to the 
WASP-10 light curve divided into individual seasons. A peri- 



odic signal was identified. Then, a sinusoid was fitted to the 
phased light curve to determine the amplitude. We deter- 
mined peak-to-peak amplitudes of 20.8 ± 1.3, 10.4 ±1.0 and 
8.5 ± 0.9 mmag for the 2004, 2006 and 2007 sea sons, respec- 
tively. Our values are similar to the results of ISmith et al] 
(|2009h who r eported 20.2 and 12 6 mmag for 2004 and 2006, 
respectively. iDesort et al] l|2007i ) showed that for a G2V- 
type star with vsini = 7 km s _1 and a photometric ampli- 
tude of 1 per cent (parameters close to those of WASP-10), 
a spot may produce RV variations with the semi-amplitude 
of ~70 m s . This estimation is consistent with the value 
determined by us from the RV model analysis. 

Stellar spots may affect the transit depth (e.g. 
ICzesla et alj|2009l ). We found in our data that transit depth 
variation is smaller than the level of 3a for individual tran- 
sits which is thus insignificant. If stellar spots are located 
along a planet's path projected onto a stellar disk, the shape 
of the transit may be deformed. This effect may have an 
influence on determining mid-tran sit times and result in a 
spurious periodic TTV signal (e.g. lAlonso et all [2009). Vi- 
sual inspection of photometric residuals (Fig. [TJ does not 
reveal features which could be related to stellar inhomo- 
geneities rather than to red noise. In particular, there are 
no deviations visible in light-curve residuals during ingress 
or egress. Such deformations caused by occulting spots lo- 
cated close to a limb of the stellar disk are expected to 
generate significant TTV amplitu de, even greater than a 
minute (|Miller-Ricci et al . 2008b). Furthermore, no peak 
close to the rotational period of the star was found in the 
periodogram of the O — C diagram (Fig. [3J) and /" v is 
not an alias of the rotational frequency of the host star 
l|Alonso et al.ll2009r i. 

The proposed model of the WASP-10 system (solution 
1) was found to be dynamically stable. Simulations were per- 
formed with the Systemic Console employing Bulirsch- 
Stoer integrator and covered 10 6 yr with a maximal time 
step of 10" 4 yr. The eccentricity of WASP-lOb was found 
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to be close to zero, while eccentricity of WASP-lOc did not 
exceed 0.05 with a median value of 0.024. 



6 AGE OF THE SYSTEM 

Given the relatively short rotation peri od of ~12 days of 
WASP-10, gyrochronology (|Barnesll2007h can roughly yield 
the a ge of the star: For a K5 dwarf with (B — V) = 1.15 
mag |Schmidt-Kalerlll982t h we get 270 ± 80 Myr only, sug- 
gesting that_th^_s_tar_^(aiid_the planetary system) is quite 
young. C hristian et alj (|2009l ) estimated the rotational age 
to be between 600 Myr and 1 Gyr by comparing the spin pe - 
riod of WASP-10 to stars in Hyades (jTerndrup et al]|2000lh 
However, one must no te that according to Fig. 7 and 8 in 
iTerndrup etTall (|2000h . there are also Pleiades members of 
spectral type around K5 (like WASP-10) with vsmi a round 
6 km/s or lower (as WASP-10. [Christian et~al] 120091 ). This 
method shows that Pleiades-like young age or an age inter- 
mediate between Pleiades and Hyades is possible for WASP- 
10. 

To confirm the relative young age of WASP- 
10, we investigated its kinematic properties. The 
proper motion of WASP-10 is (fj, a cos S, [is) = 
(21.4 ± 2.0, -28.9 ± 1.0) masyr" 1 i|Zacharias et all 

l2003h and its radial velocity is -11.44 ± 0.03 kms" 1 
( Christian et alj 120091 ') yielding a heliocentric velocity of 
(U, V, W) = (-0.3 ± 0.9, -17.3 ± 1.7, -8.1 ± 3.2) kms -1 . 
First, we compared the stars's spatial velocity with these 
of know n nearby moving gro ups including the Pleiades and 
Hyades dAntoia et all 120081 . e.g. their Figure 5, Table 3; 
IZhao et al.ll2009l , e.g. their Figure 2, Table 1). We found 
that WASP-10 cannot be associated with any of the moving 
groups. In fact, it lies in an area in the U-V plane where the 
star density is rather low, which may indicate that this star 
was either formed in isolation (very rare) or was ejected 
from its parent association. 

If WASP-10 was ejected recently, then it might be pos- 
sible to identify its parent cluster. As stars are ejected most 
likely soon after cluster formation (du e to a supernova event 



in a binary system, e.g. lBlaauwlll96ll or via dynamical inter- 
actions in a dense cluster, e.g. iPoveda et al.lll967h . the age 
of WASP-10 should be its kinematic age, i.e. it left its parent 
cluster some hundred million years ago. For such a long time, 
it is hardly possible to reconstruct the past trajectory of the 
system as well as the past trajectory for a potential parent 
cluster because both, the star and the cluster, have experi- 
enced a complete cycle (or even more) around the Galactic 
centre. Considering a time span up to 1 Gyr, even periodic 
appearances of close encounters between the star and a clus- 
ter could occur. Moreover, it is not unlikely that WASP-10 
experienced some perturbations of its path. For those rea- 
sons, the parent cluster of WASP-10 might remain uniden- 
tified if the star was ejected relatively soon after formation. 
If, for some reason, WASP-10 was ejected recently, finding 
its parent cluster could be possible. The system's peculiar 
spatial velocity is 10 ± 2 kms -1 (the heliocentric velocity 
was corrected for So lar motion and Galac tic rotation ap- 
plying the LSR from iTetzlaff et al] (|2010bl 'l and a value of 
VQ, ro t = 225 kms -1 ). Considering a typical cluster radius of 
a few parsecs, the parent cluster of WASP-10 is well iden- 
tifiable if the star system was ejected not earlier than some 



Table 5. Open clusters for which close encounters with WASP-10 
were found some million years in the past. The first column gives 
the cluster designation, the second and third columns indicate 
the separation d between WASP-10 and the cluster centre and 
the radius R of the cluster. Column four states the time r of the 
encounter (ti me before present) and column five quotes the cluster 
age as given in lDias et al.l ((2010). d was derived from 3D Gaussian 
distributions that can explain the slope of the d distribution, see 
ITetzlaff et al] ^OlOal l. 



Cluster 


d (pc) 


R (pc) 


t (Myr) 


age (Myr) 


Platais 2 


0±6 


10 


13 


400 


ASCC 123 


± 12 


6 


11 


260 



million years ago. During 20 million years, the star travelled 
about 200 pc. Taking into account the current distance to 
the Sun of about 100 pc and considering also previous age 
constraints, we selected open clusters within 300 pc from 
the Sun having ages between 100 and 800 Myr (correspond- 
ing to t he ages of the Ple iades and Hyades, respectively) 
from the lDias etahl (|2010T ) database for open clusters with 
full kinematics available with the catalogue as well as the 
Hyades cluster fWEBDAFI TMermilliod fc Paunzenllioiol '). 14 
clusters in total. We calculated the past orbits of WASP-10 
and each cluster us ing the epicycle approximation (jLindbladl 
ll959l ; IWielenlll982r i. Then, we performed a Monte-Carlo sim- 
ulation varying the distance, proper motion and radial ve- 
locity within their confidence intervals to account for their 
uncertainties. The goal was to find close encounters between 
WASP-10 and any cluster in the past. Two clusters for which 
close encounters with WASP-10 were found (within three 
times the cluster radius corresponding to approximately 3cr), 
are listed in Table [5] The cluster radius, the time of the en- 
counter and the age of the clusters are given. If WASP-10 
really originated from either Platais 2 or ASCC 123, the star 
is about 260 to 400 Myr old, in agreement with our age es- 
timate of 270 ± 80 Myr. It is, of course, more likely that the 
star is ejected early in the life of a cluster. It is generally 
not possible to trace WASP-10 backwards reliably for the 
estimated age of the system. 

The young age of the host star may explain the 
large radiu s of its planet , WAS P-lOb. The non-irradiated 
models by iBaraffe et alj (|2003T ). interpolated for such a 
young 3 Mj planet, results in a radius of ~1.2 Rj. Ir- 
radiation may increase the planetary radius by ^10 per 
cent, compared to non-irradiated models (IBaraffe et al.l 
120031 ). An analogo us procedure based on irra d iated plan- 



etary models by Fortnev. Marlev fc Barnesl |2007) and 
IBaraffe. Chabrier fc Bartmanl ( 20081 ) also gives a radius 
of ~1.2 Rj for a 300 Myr old giant planet. This es- 
timate is consistent with the r a dius of WASP-lOb re - 
por ted by IChristian et al.1 (|2009h . iDittmann et alj (|2010T ) 
and iKreicova et ah ( 2010T ) . 



7 CONCLUSIONS 

Our investigation indicates that transit timing of WASP-lOb 
cannot be explained by a constant period of the exoplanet. 



http: / /www. univie.ac.at/webda/webda. html 
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The distribution of data points in the O — C diagram reveals 
the existence of a periodic signal. We constructed a provi- 
sional hypothesis assuming the presence of the third body 
in the system. Combining results of 3-body simulations and 
RV measurements the most likely explanation of the obser- 
vations is given by a model with a second planet of mass of 
~0.1 Mj and orbital period of ~5.23 d. 

The reanalysis of RV data shows that WASP- 10b orbits 
its active host star in a circular orbit. The signature of stel- 
lar rotation was found in RV measurements. This effect was 
misinterpreted in previous studies as a non-zero eccentricity 
of the transiting planet. The second planet would generate 
wobbling of the host star with the semi-amplitude of ~14 
m s _1 - much less than the RV variation caused by the 
stellar rotation. Although we managed to distill this weak 
signal, we emphasis that further simultaneous spectroscopic 
and photometric observations are needed to study the ac- 
tivity of WASP-10 and verify the existence of the planet 
WASP-lOc. 

WASP-10 was found to be a relatively young star. As- 
suming that this finding, which is based on the stellar gy- 
rochronology, is correct, then the observed radius of WASP- 
10b is consistent with theoretical models and no tidal heat- 
ing is required to match its radius. If planetary gas giants 
form by gravitational contraction, then young planets should 
be larger than older ones. If the planets of WASP-10 formed 
at larger separation than their present location, then mi- 
gration must have taken place within the young age of the 
star. 

The existence of WASP-lOc could be independently 
confirmed if it is a transiter. Assuming that its radius is 
~0.4 Rj, i.e. similar to exopla nets of the si milar mass like 
HAT -P-llb l|Bakos et alj|20ich or Kepler-4b (|Borucki et all 
l2010l ). the expected flux drop during a transit would be ~0.3 
per cent or even greater if we consider the young age of the 
system. This could be observable with a large ground-based 
telescope. 

In the proposed model the ratio of orbital periods of the 
planets P c /Pb = 1-69 is very close to 5:3 orbital resonance. 
A significant fraction of known multi-planet systems are in 
MMRs, mainly in low-order ones. Two planets in wide orbits 
and trapped in a 3:2 MMR were found a round HD 45364 
by RV measurements l|Correia et a l. 2009). Planets in the 



4-planetary system around GJ 581 ( Mayor et alj|2009[ ) are 
very clo se to 5:2 (GJ 581b and c) and 5:3 (GJ 581e and b) 
MMRs (|Papaloizou fe TerquemlboiOl ). Therefore, we con- 
clude that the architecture of the WASP-10 system would 
not be completely unusual. 
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